Population-based study of the reproductive risk factors for transvaginal ultrasound diagnosed uterine fibroids in Nigerian women

There has been no previous systematic, epidemiological study of the reproductive risk factors for uterine fibroids (UF) in African populations despite African women having the highest burden of UF in the world. Improved knowledge of the associations between UF and reproductive factors would contribute to better understanding of the etiology of UF and may suggest novel opportunities for prevention and therapeutic interventions. We used nurse administered questionnaires to survey the demographic and reproductive risk factors of UF among 484 women who are members of the African Collaborative Center for Microbiome and Genomics Research (ACCME) Study Cohort in central Nigeria, and who had transvaginal ultrasound diagnosis (TVUS). We used logistic regression models to the evaluate associations between reproductive risk factors and UF, adjusted for significant covariates. In our multivariable logistic regression models, we found inverse associations with number of children (OR = 0.83, 95%CI = 0.74–0.93, p-value = 0.002), parity (OR = 0.41, 95%CI = 0.24–0.73, p-value = 0.002), history of any type of abortion (OR = 0.53, 95%CI = 0.35–0.82, p-value = 0.004), duration of use of Depot Medroxyprogesterone Acetate (DMPA) (p-value for trend = 0.02), menopausal status (OR = 0.48, 95%CI = 0.27–0.84, p-value = 0.01), and a non-linear positive association with age (OR = 1.04, 95%CI = 1.01–1.07, p-value = 0.003). Other reproductive risk factors that have been reported in other populations (age at menarche and menopause, and oral contraceptives) were not associated with UF in this study. Our study confirms some of the reproductive risk factors for UF that have been found in other populations and shows that some of them are stronger in the Nigerian population. The associations we found with DMPA suggest opportunities for further research to understand the mechanisms of action of progesterone and its analogues in the etiology of UF, their potential use for prevention and treatment of UF.

of UF, the association with oral contraceptives remains unsettled.Some papers report inverse associations while others found positive or no associations 1,7,9,15 .Hormone replacement therapy in postmenopausal women was associated with increased risk of UF in some studies 7 .Studies have also shown interactions between reproductive risk factors and race in admixed societies like the US 7,16 .While the growth rate of UF tends to decline in white women as they approach menopause, similar decline has not been observed in black women 16 .
In this paper, we report our adequately powered, systematic, population-based study of the associations between reproductive risk factors and transvaginal ultrasound diagnosed UF in Nigerian women.

Methods
The ACCME cohort study profile and the characteristics of the study participants have been previously described in detail 17 .We enrolled 11,700 women who were at least 18 years old at baseline, HIV negative, and had a history of penetrative vaginal intercourse but no previous history of cervical abnormalities, cervical cancer, or hysterectomy into the ACCME cohort from 2014 to 2016.Trained research nurses administered validated questionnaires to collect demographic, lifestyle, obstetric and gynecologic history, medical history, sexual behavior and practices, diet, and physical activities.
In 2020, for this study, we developed and internally validated a supplemental questionnaire among 20 ACCME participants, who were not included in the rest of the UF study.The questions were: a. Has a doctor ever told you that you have uterine fibroid and you have not received any treatment for it?b.Has a doctor ever told you that you do not have uterine fibroid?c.Have you ever been told that you have a uterine fibroid on abdominal ultrasound scan and have not received treatment for it?d.Have you ever been told that you do not have uterine fibroids on abdominal ultrasound scan?e. Do you have any symptom-excessive bleeding during menstrual periods, severe pelvic/abdominal pain, inability to conceive, abdominal swelling-that have been attributed to uterine fibroids?
Using random numbers generated with Stata software, we identified and invited 486 women from the ACCME cohort to participate in this study.All those invited responded and consented.We retrieved the participants' epidemiological data from the ACCME Cohort baseline dataset.The variables we selected included age in years, years of education completed (no formal education, 1-6 years of schooling, 7-12 years of schooling, any schooling after 12 years including university education, higher diplomas, etc.), occupation, religion, and marital status (married, unmarried).To compute socio ® -economic status (SES), we calculated the 'wealth index' using the following variables: house ownership and type of house owned (e.g., home, apartment, house, or duplex); source of drinking water (e.g., from outside, well, borehole, piped or bottled); type of cooking fuel; use of separate room for cooking; type of toilet; and ownership of household goods including car and refrigerator.We used principal component analysis (PCA) with varimax rotation to compute factor scores based on the sum of responses to these variables weighted by their factor loading.We used the first component in the PCA that explained 35% of the variation in the data to generate a wealth index (39).The wealth index variable was used to classify participants into low SES (lowest 50% of the score distribution), middle (middle 30%) and high (highest 20%) socioeconomic classes.
We also retrieved data on reproductive risk factors including age at menarche, age at sexual debut, parity, history of spontaneous or induced abortion, history of ectopic pregnancies, age at menopause, menopausal status (premenopausal vs. postmenopausal-individuals who reported that they were not sure of their menopausal status were classified as post-menopausal in multivariable analyses), stillbirths, history of regular contraceptive use (yes vs. no), age at commencement of regular contraceptive use, and duration of use of either oral contraceptives, Depot Medroxyprogesterone Acetate (DMPA), or contraceptive implants.
Participants were offered transvaginal ultrasound scan (TVUS) performed by clinical radiologists whose fellowship was certified by the West African College of Surgeons.UF identification was based on the Muram criteria and included all relatively spherical, echogenically distinct masses in the myometrium that were more than 0.5 cm in diameter 18 .Participants with UF on TVUS were referred to the gynecologic clinic for standard of care treatment.
We used the Table 1 command in Stata 18 ® to generate our table of descriptive statistics showing distribution of variables among cases and controls, and p-values for chi-square tests of categorical variables and t-tests of continuous variables comparing women with UF and those without.To construct our regression models, we identified variables with p-values < 0.20 in age-adjusted analyses and included them in our multivariable logistic models.We used restricted cubic splines to explore continuous variables such as age at menarche, age at sexual debut, and age at menopause for potential non-linear relationships with risk of UF.In the multivariable models, we used Wald test to identify covariates with significant associations with risk of UF at p-value < 0.05.We analyzed all data using Stata 18 ® (College Station, Texas 77845 USA).

Patient and public involvement
Patients and the public were not involved in the design and analysis of this project.The study is population based and thus did not enroll individuals who presented as patients in a medical institution.The findings on TVUS were shared with participants who were counselled on accessing the standard of care for UF in Nigeria.The idea and design of the study was based on clinical experience with patients, but these patients were not involved in the recruitment or conduct of this study.The ACCME Project has an active website and social media that is used to engage study participants and disseminate the results of this study.This website is also used to grant controlled access to data used in this study.

Results
Of the 486 women invited, there was complete data for analysis in 484 women.The mean (SD) age of all participants in this study was 37.0 (9.20) years, and this was not significantly different comparing women with UF (mean = 36.9,SD = 9.03) with those without UF (mean = 37.1, SD = 9.35) (p-value = 0.88).Women with UF were more likely to have received more formal schooling (p-value = 0.04) and have single marital status (p-value = 0.048) compared to women without UF (Table 1).

Reproductive risk factors and UF (Table 2)
Women without UF had more children (mean = 3.12, SD = 1.97 compared to mean = 2.62, SD = 1.99; p-value = 0.006), were less likely to be nulliparous (p-value = 0.004), more likely to have had any type of abortion (p-value = 0.02), more likely to have regularly used any form of contraception (p-value = 0.02), and were more likely to have used DMPA for at least a year or more (p-value < 0.001).There were no significant associations between age at menarche, age at sexual debut, age at regular use of any form of contraception, history of ectopic pregnancy, age at menopause, menopausal status, stillbirths, duration of use of oral contraceptives, and duration of use of contraceptive implants, and TVUS diagnosed UF.

Multivariable analyses (Table 3)
In age-adjusted logistic regression models, parity, menopausal status, history of any type of abortion, history of regular use of any form of contraception, and duration of use of DMPA were significantly associated with www.nature.com/scientificreports/confounding.We conducted likelihood ratio test of the model with SES and that without SES and there was no significant difference (p-value = 0.73), so we removed SES from our final model.Age became strongly associated with UF in multivariable logistic regression models whereas it was not associated when it was regressed against UF thus suggesting that the association may be nonlinear, or the association with age became more obvious after adjustment for confounders.To check the linearity of the relationship between age and UF, we estimated a locally weighted scatterplot smoothing (LOWESS) curve of the relationship between age and UF in the logistic regression model 19 .Fig. 1 shows the graphical relationship suggesting increasing association with UF till age 35 followed by a shallow decline and a subsequent rise in risk of UF with increasing age.

Discussion
In this first systematic epidemiological study of the associations between reproductive risk factors and TVUS diagnosed UF in a population-based cohort of Nigerian women, we found associations with parity, abortion, regular use of DMPA for contraception, menopausal status, and a non-linear association with age.Other reproductive risk factors that have been reported in other populations (age at menarche and menopause, and oral contraceptives) were not associated with UF in this study 1,2,7,13 .Previous studies of reproductive risk factors for UF in African populations showed associations with age and parity [20][21][22] .However, none of the previous studies in Africa were population-based and they did not use validated tools for diagnosis of UF.They were also often poorly designed and inadequately powered 2 .
Several studies have shown an association between parity and UF but this is the first study to show similar and stronger association in African women compared to women with non-African ancestries 7,23 .In this study, we found 17% reduction in risk of UF for each child and 59% reduction in risk of UF comparing parous Table 3. Multivariable analyses of reproductive risk factors for TVUS diagnosed UF, Nigeria, 2021.a These variables were replaced for each other in the multivariable models and both were not simultaneous included in same models.b Variables did not achieve significance in multivariable models and were dropped.◊ p-value for trend.www.nature.com/scientificreports/with nulliparous women 7 .These risk reductions are substantially higher than the 20% to 50% risk reductions previously reported from studies conducted in other populations 7 .There are multiple potential mechanisms whereby parity may influence the development and growth of UF.These include direct effect of pregnancy associated hormones on the myometrium and the extracellular matrix that constitutes UF, and non-hormonal mechanisms 7,10,11,15,24 .Fibroids may also involute during the remodeling of the uterus, post-partum 25 .

Age
The history of any type of abortion was associated with reduced risk of UF by about 17% in this study.This effect may be independent of the effect of parity.Previous studies have reported conflicting associations between abortions and UF 7,24,26 .While some have suggested that a causal relationship between spontaneous abortions and UF is questionable, there is little doubt that ascertaining causality or reverse causality in this association would be challenging, even for the most robustly designed and implemented study 27 .This finding may be a reflection of the association between pregnancies and UF in general and not specific effect of abortions on risk of UF.
We found a marginal association between regular use of any type of contraception and UF, but this did not persist in multivariable models.This may be because more granular data on type and duration of use of different contraceptives than what is available in our study is required to adequately capture this exposure.There was no association with use of oral contraceptives.Most of the previous studies of oral contraceptives and UF have shown either decreased or null associations 7,28,29 .There are serious methodological challenges with research on use of oral contraceptives in sub-Saharan Africa where uptake is low, there are variations in the compositions of oral contraceptives sold, and adherence with regular intakes is uncertain 30,31 .We found significant association with use of DMPA and this declined with increasing duration of use in multivariable models.This is similar to findings in other populations 7,15,[32][33][34] .DMPA may be associated with reduced risk of UF due to estradiol mediated reduction in progesterone receptor expression or differences in the mechanisms of action of exogenous compared with endogenous progesterone 34 .
The risk of UF in our study rose rapidly from early adulthood to about 35 years of age after which it declined slightly, then it resumed a gentler increase in risk.This finding is consistent with the pattern of change in reproductive risk factors with aging and the findings from other studies, including cohort studies that reported peak age incidence of 40-44 years 35 .The earlier peak age in this study may be a reflection of the population distribution of the underlying study population and is not an age incidence, therefore, it should not be compared with results from cohort studies.Nevertheless, recent global studies of UF suggest that the associations of UF with age peaked at 35 to 44 years, then declined with advancing age 36,37 .There was a 53% reduction in the risk of UF at menopause in this study.The natural regression of fibroids at menopause has been established in several studies and this has been attributed to the changes in estrogen and progesterone levels that occur at menopause 7,38 .
Despite our large sample size, population-based enrollment of study participants, and use of TVUS, our study has some limitations.We used data from questionnaires that were administered before we enrolled participants for the UF study.While this may reduce recall bias, the life course of UF is dynamic and may have changed in individual participants between enrollment into the ACCME study and TVUS.Large, prospective cohort studies in African populations are required to adequately account for secular trends in life-course of UF and their relationships with risk factors.We did not have data on all the reproductive risk factors for UF.Future studies focused on UF may seek to elicit more detailed responses but at the risk of over-burdening participants and reducing the quality of responses.

Conclusions
In this study, we confirmed many of the associations between UF and reproductive risk factors that have been described in other populations for the first time in a population-based, study of TVUS diagnosed UF in an African population.This is important because African women bear the highest burden of UF in the world.Studies like ours improve understanding of the risk of UF which may lead to identification of biomarkers and therapeutic targets.They also create opportunities for early detection and preventive interventions.

Table 1 .
Baseline characteristics of women with TVUS diagnosed UF and those without TVUS diagnosed UF, Nigeria, 2021.
Vol:.(1234567890) Scientific Reports | (2023) 13:18926 | https://doi.org/10.1038/s41598-023-44703-5www.nature.com/scientificreports/history of TVUS diagnosed UF.However, only age, history of any type of abortion, duration of use of DMPA, and menopausal status were significantly associated with UF in multivariable logistic regression models.Marital status and history of regular use of contraceptives did not retain the associations observed in age-adjusted models.Education was marginally associated with UF in age-adjusted models but not in multivariable models.While socio-economic status (SES) was not significantly associated with UF in age-adjusted and multivariable models, we retained it in the multivariable models because of its potential to account for unmeasured